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ABSTRACT 

We report on a Spttzer/IRS (mid-infrared), RXTE /PC A+REXTE (X- 
ray), and Ryle (radio) simultaneous multiwavelength study of the microquasar 
Cygnus X-1, which aimed at an investigation of the origin of its mid- infrared 
emission. Compact jets were present in two out of three observations, and we 
show that they strongly contribute to the mid-infrared continuum. During the 
first observation, we detect the spectral break — where the transition from the 
optically thick to the optically thin regime takes place — at about 2.9 x 10^^ Hz. 
We then show that the jet's optically thin synchrotron emission accounts for 
the Cygnus X-l's emission beyond 400 keV, although it cannot alone explain 
its 3 — 200 keV continuum. A compact jet was also present during the second 
observation, but we do not detect the break, since it has hkely shifted to higher 
frequencies. In contrast, the compact jet was absent during the last observation, 
and we show that the 5 — 30 /im mid- infrared continuum of Cygnus X-1 stems 
from the blue supergiant companion star HD 226868. Indeed, the emission can 
then be understood as the combination of the photospheric Raleigh- Jeans tail and 
the bremsstrahlung from the expanding stellar wind. Moreover, the stellar wind 
is found to be clumpy, with a filling factor /oo ~ 0.09 — 0.10. Its bremsstrahlung 
emission is likely anti-correlated to the soft X-ray emission, suggesting an anti- 
correlation between the mass-loss and mass- accretion rates. Nevertheless, we do 
not detect any mid-infared spectroscopic evidence of interaction between the jets 
and the Cygnus X-l's environment and/or companion star's stellar wind. 

Subject headings: binaries: close — X-rays: binaries — Infrared: stars — dust, 
extinction — Stars: individual: Cygnus X-1 — ISM: jets and outflows 
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Introduction 



Accretion— ejection phenomena are respo nsible for the strong 



correlations observed in microquasar s (see e.g. 



Gallo et al 



2003 



Coriat et al 



Mirabel et al 



1998 



radio / infrared /X-rays 



Corbel fc Fender 



2002 



20091 ). Indeed, some of the accreted material is returned 
to the interstellar medium (ISM) through collimated relativistic ejections, whose power 
law-like synchrotron radiation {F^, oc i/") is thought to extend from the radio to the X-ray 
domains. The best known are the so-called compact jets, only detected in the hard state 
(HS), whose spectral continuum is well-modeled by a fiat or weakly- inverted power law in 
the optically thick regime (0 < a < 0.2) and decreasing in the optically thin one, with 
a ranging between — 0.4 and —1 depending on the electron energy distrib ution (for a 



-0.6 canonical value, 



Blandford fc Koenigl 



1979 



Falcke fc Biermann 



1995|). To date, the 



spectral break where this transition takes pla ce has only been detected photometrically 



in two sources, GX 339-4 and 4U 0614+091 fICorbel fc Fender 



2002 



Migliari et al. 



2006 



20101). Yet the knowledge of this cut-off frequency is crucial for a full understanding of the 
properties of the compact jets, as it is closely related to the original physical conditions, 
such as the magnetic field, the base radius of the jet, and the total energy of the electron 
population. 



Among microquasars, Cygnus X-1 (IBowyer et al. 



19651 ) is one of the few classified 



as a high mass X-ray binary (HMXB). It is composed of an 8—16 Mq black hole 



(jShaposhnikov fc Titarchuk 



2007; 



09.7 ab supergiant HD 226868 flWalborn 



Caballero-Nieves et al 



It has a quasi-circular 5.6 d orbit ( iBrocksopp et al 



1973 



lorio 



2008 



1999a 



2009) bound to the 17-31 



Caballero-Nieves et al 



20091) 



, and re : 



assoc iated inclination in the range of 23° — 43° fiGies et al. 



2003 



erence therein), with a n 



Caballero-Nieves et al. 



20091 ). The distance to the sour ce is still a matter of debate, but severa. 



values between 1.5 and 2.5 kpc (jZiolkowski 



2005 



Caballero-Nieves et al. 



studies have found 



20091. and reference 
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therein). Recently, IXiang et al.l (l201l|) found 1 .80 kp c with a 7% accuracy from X-ray dust 



scattering halo studies. Moreover, iReid et al 



(120 111 ) derived 1.86 kpc with a 6% accuracy. 



from radio parallax measurements obtained with the very large baseline array (VLBA) at 
8.4 GHz between 2009 Jan 23 and 2010 Jan 25. Both results are consistent with 2 kpc 
within the uncertainties, and we adopt this value in this study. 



Cygnus X-1 spends two-thirds of its time in the HS (jWilms et al. 



20061). Indeed, 



although the companion is close to filling its Roche lobe, the transfer of mate rial takes place 



through the strong ste 


liar winds that are 


1986; 


Gies et al. 


2003; 


Hanke et al. 


2009) 



Gies fc Bolton 



20091 ). This lower accretion rate transfer is less effective 



than for Roche lobe overflow microquasars, there by most of the time preventing the tran 



Pottschmidt et al. 



2000 



2003h . Instead, 



sition to a disk-dominated soft state (SS, see e.g. 
the source often exhibits a so-called failed-state transition, characterized by a continuous 
increase of the X-ray flux up to an intermediate state (IS), followed by a strong decrease 
and a return to the HS, with possible launches of discrete ejections. Finally, its emission 
appears to be in phase with t he orbital period in most of the spectral domains in which 



the source has been detected, ( Pooley et al. 



1999 



Brocksopp et al. 



1999a 



Lachowicz et al. 



20061 ). This includes the radio activity, which is similar to that of other microquasars. 



In particular, the compact jets, characterized by a very fiat spectrurn (a < 0.1) and a 



10—25 mJy flux density level are detected in the HS (iFender et al. 



2000 



Stirling et al 



200l|). 



In this paper, we report on a mid- infrared (MIR), radio, and X-ray multiwavelength 
study of Cygnus X-1, with a focus on MIR spectroscopy. Our goals are to (1) assess the 
origin of the MIR emission, (2) compare the possible variation of the MIR continuum to 
that of the X-ray and/or radio emission, and (3) detect any spectroscopic feature arising 
from an interaction between the jets, the stellar wind, and/or the interstellar medium 
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(ISM). The observations and the data reduction are presented in Sect. 2, while Sect. 3 
is devoted to the X-ray and radio study. In Sect. 4 and Sect. 5, we analyze the broad 
band MIR to radio SEDs of the source and the spectroscopic contents of the MIR spectra, 
respectively. Finally, we discuss the outcomes in Sect. 6. A companion paper, Lee et al. (in 
prep), will focus on the dust content of Cygnus X-1. 



2. Observations and data reduction 



Cygnus X-1 was observed on UT 2005 May 23, June 7, and July 2 (P.I. S. Heinz) with 



the Infra Red Spectrograph 



{Spitzer; 



Werner et al. 



(IRS; 



Houck et al. 



2004 ) on-board the Spitzer Space Telescope 



2004J ). In the X-ray domain, the source was observed simultaneously 



with the Proportional Counter Array (PCA; 



Jahoda et al 



Timing Experiment (HEXTE; iRothschild et al 



20051 ) and the High Energy X-ray 



19981 ) instruments (P.I. J. Wilms), and we 



made use of 



:he 1. 2—12 keV public data provided by the All-Sky Monitor instrument (ASM; 



Levine et al. 



19961 ). all mounted on the Rossi X-ray Timing Explorer (RXTE). Finally, 
Cygnus X-1 was regularly monitored at 15 GHz with the Ryle radio telescope (P.I. G. G. 
Pooley) and we used the 20 s and 5 min light curves that are simultaneous with the IRS 



spectra, for the analysis presented in this paper (see Tabled]). 



2.1. IRS observations 

We used the SL2 (5.20 - 7.70 ^um), SLl (7.40 - 14.50 /xm), LL2 (14.00 - 21.30 /im), and 
LLl (19.50 — 38.00 fim) modules at two nodding positions for sky subtraction. The total 
exposure times were each time set to 60 s in SLl and SL2, 120 s in LL2, and 2x120 s in 
LLl. Basic Calibration Data (BCD) were reduced following the standard procedure given 
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in the IRS Data Handbooljl]. The basic steps include bad pixel correction with IRSCLEAN 
vl.9, sky subtraction, as well as extraction and calibration (wavelength and flux) of the 
spectra — with the Spectroscopic Modeling Analysis and Reduction Tool software SMART 
v8.1.2 — for each nodding position. Note that in LL2 and LLl, where the sky appeared 
weakly non-uniform between both nods, the subtraction residuals were fitted by a first order 
polynomial. Spectra were then nod-averaged to improve the signal-to- noise ratio (SNR). 

The flux calibration appeared consistent between the different modules at Obs. 1 and 
Obs. 2, with differences less than 5%, and this allowed us to combine all the spectra by 
scaling them to SLl, which has the highest flux level with respect to SL2 and LL2. However, 
at Obs. 3, SL2 had a flux level about 14% lower than that of SLl, while no such problem 
was detected with LL2 nor LLl. We were not able to definitively assess the reason for such 
a discrepancy, but a misalignement of the SL2 slit is the most reasonable explanation. We 
therefore scaled the SL2 flux level to that of SLl. 



2.2. RXTE observations 



In the reduction of the data from the two pointed instrumen ts PGA and HEXT E, we 



Wilms et al. 



20061). Data 



follow the procedure used in previous studies of Cygnus X-1 (e.g. 
were extracted using HEASOFT 6.9 for all times at least 10 minutes after passages of the 
South Atlantic Anomaly. During the observations, two (first two RXTE observations) and 



three (third RXTE pointing) proportional counter units were in operation. 



http:/ /ssc. spitzer.caltech.edu/irs/irsinstrumenthandbook/IRS _Instrument_Handbook.pdf 
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2.3. MIR data dereddening 



Our Cygnus X-1 observations show a strong absorption feature around 9.7 fim. At 
the distance to the source (~ 2 kpc), we expect a considerable column of interstellar 
material, so this feature is most likely from silicate dust in the ISM along the line-of-sight. 
Several consistent measurements of the visible extinction are found in the literature (see 

and reference therein). In the companion paper (Lee et al., 



Caballero-Nieves et al. 



2009 



in prep), we derive an absolute visual extinction Ay ^ 2.95 from t 



silicate absorptio n featu re rg.y, defined by Ay = 18.5 x rg y ( jPraine 



that of 



Wu et al 



le op acity of the ISM 



20031 ). consistent with 



fll982[ ). Ay = 2.95 ± 0.22 with Ry = 3.1. 



We therefor e use t his value to deredden the spectra with the extinction laws given in 



Chiar fc Tielend ( l2006l ). one of the most recent works covering the entire IRS spectral range. 
In their paper, the authors derived the A\/Ay^ ratio rather than the usual A\/Ay one. To 
express the extinction ratio in the standard way, we assign to the value derived from the 



law given in iFitzpatrickl fll999f ) for the diffuse ISM {Ry = 3.1), which is Ak = 0.111 x Ay. 
Moreover, these authors find that beyond 8 fim, where silicate absorption dominates, the 
MIR extinction is different for the diffuse ISM and the Galactic center. They consequently 
propose two extinction laws, and we explored both options for dereddening our spectra 
of Cygnus X-1. We find that the Galactic center case clearly overestimates the silicate 
absorption depth, and we therefore adopt the corrections based on the extinction due to the 
diffuse ISM. We note that a single extinction correction (i.e.. Ay = 2.95 ± 0.22) produces 
equally good results for all three epochs. This gives us further confidence that the observed 
silicate absorption features are caused by dust in the diffuse ISM as opposed to material 
directly associated with the Cygnus X-1 system. 

The resulting dereddened spectra displayed in Fig. [1] show unambiguously that the 
MIR emission of Cygnus X-1 is variable. Indeed, the MIR continuum was brighter during 
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Obs. 1 and Obs. 2, the difference with respect to Obs. 3 increasing with wavelength, from 
no difference at 6 /^m to an excess of 15—20 mJy beyond 10 /xm. 



3. X-ray/radio behavior 

Since dust heating as a contribution to the MIR emission can be ruled out due to 
the lack of observed characteristic dust spectral features (Lee et a., in prep), we consider 
four main alternative physical processes: (1) the stellar photosphere, whose emission is 
constant, (2) bremsstrahlung from the expanding stellar wind, (3) synchrotron emission 
from jets, (4) bremsstrahlung from the accretion disk wind. We use the simultaneous radio 
and X-ray data to assess the extent to which the last two components are present during 
the observations (see Sect. 4 for details). 



3.1. RXTE spectra 



The RXTE spectra are fitted between 3 and 200 ke V in XSPEC v. 12.5.1. We tested 



several continuum models s uch as the stand ard a— disk (IMitsuda et al. 



Comptonization component (ITitarchuk 



19841 ) plus a 



19941 ) modified by photoelectric absorption, but we 



find that all the spectra are best described by the phenomenological model combining a 
broken power law with a high energy cut-off and a gaussian to account for the iron emission, 
both modified b y photoelectric ab sorption. These results are consistent with previous 



studies (see e.g. 



Wilms et al. 



20061 ). The best-fit parameters for all the observations are 



listed in Table [2] and the fits are displayed in Fig. El 



Based on the classification of Cygnus X-l's spectral states as given by 



Wilms et al. 



( l2006l ). these fit results show unambiguously that the source is in the HS during Obs. 1 



and Obs. 3, with Fi photon indices of about 1.83 and 1.75, respectively. The 3 — 200 keV 
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unabsorbed flux is nevertheless about 50% larger during Obs. 1 than during Obs. 3 which, 
considering the X-ray/radio correlation observed in microquasars in the HS, suggests a 
stronger contribution of compact jets to the overall emission of Cygnus X-1 during Obs. 1. 
In contrast, the spectrum of the source during Obs. 2 is softer, with Fi ^ 2.08. Such a 
photon index is consistent with an IS. But similarly, the high 3 — 200 keV unabsorbed flux 
suggests a stronger contribution of compact jets than during Obs. 3. 



3.2. Light curves 

In Fig. 131 we display the Cygnus X-1 radio (15 GHz, 5 min time resolution) and ASM 
(1.2—12 keV, day-averaged) light curves, as well as the ASM HR2 (C/A) hardness ratio, 
between MJD 53500 and MJD 53580, bracketing our observations. Clearly, it is seen that 
the count rate is lower and the hardness ratio higher during Obs. 3, i.e. when the MIR 
continuum appears to be at its lowest level (see Fig. [H blue spectrum). Consistent with the 
spectral fitting, the ASM light curve and HR2 hardness ratio show that the source during 
Obs. 2 is in a softer state and exhibits a slightly higher fiux than during Obs. 1 (see also 
Fig. Hj). This is in good agreement with the light curve trend shown in Fig. |3] between 
MJD 53510 and MJD 53530, which shows a sharp decrease in the hardness ratio and a 
slow increase in the count rate. This probably points to a higher accretion disk activity 
becoming unsteady around MJD 53520. Then, the count rate drops between MJD 53530 
and MJD 53550 before stabilizing, while the hardness ratio follows the oppo site path. 



his b ehavior is characteristics of a failed-state transition as characterized by 



fcooeh . 



Wilms et al. 



Moreover, the Ryle light curves during Obs 1. and Obs. 2 exhibit an average radio flux 
between 12 and 25 mJy, while the radio flux is on average lower and less steady during 
Obs. 3. This behavior is shown by the 20 s time resolution Ryle light curves, displayed in 
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Fig. ini which cover the exact time range of each IRS observation. It is clear that Obs. 1 
and Obs. 2 are consistent with compact jets. In contrast the radio activity during Obs. 3 is 
very unstable, covering a range from non-detection to 48 mJy. This behavior is not by itself 
inconsistent with the presence of compact jets, and possible explanations of this flickering 
could be attributed to some instabilities internal to the jet or some peculiar accretion flow 
properties. 

Finally, it is worth mentioning t hat there have been indi c ations, from the tim ing 



behavior, of two different HS regimes (jPottschmidt et al. 



2003 



Axelsson et al 



20081 ). One 



of them is characterized by a very high hardness level, a low X-ray flux, and characteristic 
frequencies shifted to the lower edge in the power density spectrum (PDS). This is clearly 
the case of the HS detected during Obs. 3 (see Grinberg et al. in prep, for the full PDS), 
which is actually th e beginning of the hardest X-ray activity ever observed for Cygnus X-1 



( iNowak et al 



201l[ ). In contrast, the HS properties during Obs. 1 are consistent with what 



is expected in the canonical picture. 



4. The MIR continuum fitting 

We use spectral fltting to understand the origin of the continuum variation between 
the three IRS observations. Between nod 1 and nod 2, we on average found 4% and 5% 
flux differences in SL and LL modules, respectively. These quantities were quadratically 
added as systematic errors to take into account uncertainties due to flux calibration. Since 
Cygnus X-1 was never in the SS during our observations, we rule out bremsstrahlung from 
the accretion disk wind, and we only consider thermal emission from the companion star's 
photosphere, bremsstrahlung from the expanding stellar wind, and/or synchrotron from the 
compact jets. Moreover a bump is present in the three spectra between about 17 and 22 fim. 
The origin of this bump is a matter of debate, and it could be due to the dereddening 
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process or to irradiation of dust around the BH 
more likely to the Spitzer /IRS LLl 24 deficit 



S. Markoff, private communication), but 



4.1. Obs. 3: the stellar continuum 



Considering the very unstable radio activity of Cygnus X-1 during Obs. 3, a jet 
contribution to the MIR spectrum is very unlikely. We therefore only consider the emission 
from the star. To model the photosphere, we simply use a black body. Cygnus X-l's 
companion star, HD 226868, is an 09. Tab supergiant and its unabsorbed emission peaks 
in the ultraviolet (UV). The MIR alone is therefore not sufficient to constrain the star's 
temperature, so we fix it to the most recent derived value, i.e. 28000 K as given in 



Caballero-Nieves et al 



(120091) ■ based on theoretical 



their seminal papers, I Wright fc Barlow! (Il975[ ) and 



fitting to the u l travio let spectrum. In 



Panagia fc Fellil (119751 ) showed that the 



bremsstrahlung from a homogenous spherically expanding wind would be responsible for a 
MIR and radio excess oc z/°'^ in supergiant stars. We therefore add a power law to take the 
wind emission into account. This simple model is summarized by Eq. [1] 



fM=(|)%B(..28000) + iV,(^)"' (1) 

where i?^, and are the companion star's radius and distance, B is the black body 
emission, and and are the bremsstrahlung fiux density at 15 GHz and spectral index. 

The best fit is displayed in Fig. E] and the best-fit parameters are listed in Table [31 
The ratio R^/D^ is consistent with the one expected from the Cygnus X-l's companion 
star, since the derived stellar radius for a typical 2 kpc distance is about 19.9 i?©, the 
typical value for a typical late O supergiant. Moreover, the power law spectral index is 



^http: / / ssc.spitzer.caltech.edu/irs/features/ #3_LLl_24u_Deficit 
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characteristics of bremsstrahlung from an expanding stellar wind, with at ~ 0.55, consistent 
with the 0.6 theoretical value within the 90% uncertainties. The MIR continuum of 
Cygnus X-1 during Obs. 3 can therefore be described by that of HD 226868. This alone 
does not rule out the possible presence of compact jets. Nevertheless, both the average 
X-ray and radio fluxes are at least 50% lower than during Obs. 1 and 2. So, even if they are 
present, either the compact jets are too faint to be detected in the MIR, or they break to 
the optically thin regime beyond 30 fim. Indeed, in this reginie, the IR flux is expected to 



vary with respect to the radio flux as Fjr ^ Fradio^^^^'^ (Heinz 



20041 ). Considering the low 



radio flux level, this very probably excludes a significant contribution of a compact jet to 
the MIR spectrum of Cygnus X-1 during Obs. 3. 



4.2. Obs. 1 and Obs. 2: MIR contribution from the compact jets 

The radio behavior does exhibit compact jets that may be responsible for the observed 
continuum increase. We consider two cases to model their contribution, (1) the spectral 
break does not occur in the IRS spectral range so the optically thick synchrotron emission 
can be described by a simple power law, or (2) the spectral break is located in the IRS 
spectral range and the emission of the compact jets can be modeled by a broken power law. 

To better constrain the compact jets emission parameters, we include the average 
15 GHz flux densities of Cygnus X-1 during the IRS/LL observations. We then fit these 
MIR/radio spectral energy distributions (SEDs) by fixing the black body parameters and 
bremsstrahlung spectral index to those of Obs. 3. 



4-2.1. Case 1: simple power law 



The overall equation describing the Cygnus X-1 emission is: 



- 13 - 



f (.) . (|)% BW, 28000) + A', {^y- + N, (2) 

where Nj and ctj are the synchrotron flux density at 15 GHz and spectral index, respectively. 

Here, was left free to mimic a possible change in the bremsstrahlung intensity due 
to a variation of the star's mass loss rate. The best-fit parameters are listed in Table HI 
and the corresponding fitted SEDs displayed in Fig. [71 Both fits are consistent with a 
contribution of compact jets to the MIR emission of Cygnus X-1 during Obs. 1 and Obs. 2. 
Indeed, the derived spectral indices are 0.01 ± 0.02 and 0.05 ± 0.01, respectively, consistent 
with typical values. Moreover, the bremsstrahlung flux density at 15 GHz during Obs. 1 is 
found to be similar to that during Obs. 3 (the no-jet case, see Sect. 4.1), 0.66 ± 0.06 mJy 
compared to 0.74^q q7 mJy. In contrast, Ng during Obs. 2 is lower (0.44^qq5 mJy). 



4.2.2. Case 2: broken power law 
Here, the overall equation describing the Cygnus X-1 emission is: 



FM=(|)'xS(-.28000)+iV.(^)°%. 



(3) 

where N^, a^i, aj2, and are the jets' synchrotron flux density at 15 GHz, optically thick 
and thin spectral indices, and cut-off frequency, respectively. 

When we left free to vary all the parameters of the broken power law and the 
bremsstrahlung normalization Ng at the same time; however this results in poorly 
constrained fits due to several degeneracies between the various fit parameters. We therefore 
proceeded with three scenarios, each time fixing one parameter only: (a) is fixed to —0.6, 
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i.e. the theoretical spectral index expected from the optically thin synchrotron emission of 
the compact jets, (b) is fixed to the values found in case 1, (c) A^^ is fixed to the value 
found during Obs. 3, i.e. we assume that the bremsstrahlung intensity is the same for all 
the observations. The best-fit parameters are listed in Tables [5l [6l and [71 and the SEDs are 
displayed in Fig. [HI respectively. 

4-2.3. Interpretation 

A broken power law gives better results to describe the emission of the compact jets 
during Obs. 1, with reduced values of 0.81, 0.85, and 0.78 for scenarios (a), (b), and 
(c), respectively, versus 0.96 for a simple power law. The low reduced ^ire largely 
due to artificially forced high IRS systematic errors that were quadratically added to the 
uncertainties. Nevertheless, Ax^ > 34 for only one DOF less, and the three fits give very 
similar results. The spectral break is clustered in the MIR range (2.70 — 2.94) x 10^^ Hz, 
i.e. about 10.06 — 11.11 /im, and the optically thin spectral index for (b) and (c) is 
always consistent with the expected value (between —1 and —0.4) within the uncertainties. 
Moreover, for (a), the bremsstrahlung intensity is found similar to that during Obs. 3, a 
hint that the mass-loss rate is the same during both observations, as expected in the HS. 
So, based on the aforementioned reasons, we conclude that the consistency of all the fits 
strengthens the need for a broken power law to describe the MIR continuum of Cygnus X-1 
during Obs. 1. 

The interpretation of the results for Obs. 2 is more complex, even though the broken 
power law gives lower reduced (0.69 against 0.72). First, Ax^ ^ 9.7, which is far less 
significant than the values found for Obs. 1. Moreover, the fits with the broken power 
law are not physically relevant. Indeed, the fit for scenario (b) gives a positive photon 
index after the spectral break {a ~ 1.21), which is clearly inconsistent with a compact 
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jet. The two other fits are consistent with compact jets, but give cut-off frequencies that 
are about twice smaller than during Obs. 1. This is counter-intuitive since the opposite is 
expected. Indeed, the source during Obs. 2 is not in the HS but in a failed-state transition, 
i.e. an IS with a softer and brighter spectrum. In models of optically thick jet emission, 
the break frequency is coupled to the jet luminosity: if the luminosity increases, the break 
frequency also increases. An increase in the optically thick flux requires an increase in the 
jet energy density and/or the cross section of the jet, which always increases the optical 
depth, thus increasing the break frequency. The fact that our best flts for Obs. 2 indicate 
the opposite would imply a change in the jet geometry, e.g., a signiflcant increase in the 
distance of the base of the jet along the jet axis or an increase in the opening angle of the 
jet (such that the footprint of the base of the jet shrinks, decreasing its optical depth). 
Given the level of uncertainty in our data (demonstrated by the fact that our models are 
under-constrained), such conclusions are still rather speculative at this point. Moreover, 
several studies of Cygnus X-1 sho wed the existenc e of an anti-correlation betw een the 



mass-loss and the accretion rates (IGies et al. 



2003 



Boroson fc Dil Vrtilek 



2010l ). maybe 



due to photoionization of the wind by the soft X-ray emission. During the IS and SS, the 
bremsstrahlung contribution should therefore be lower than in the HS. This is actually 
exactly the result of case 1, where we consider a simple power law instead of a broken one. 
The bremsstrahlung intensity is found to be about 40% lower than that during Obs. 3, 
while the cut-off frequency is assumed to occur below 5.48 fim. We believe that it better 
describes what occurs during Obs. 2. 



5. MIR spectral lines analysis 

In the companion paper (Lee et al., in prep), we show that the presence of warm dust 
around Cygnus X-1 is very unlikely and that the strong silicate absorption feature is due 
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to the ISM. Here, we study the spectral hne variabihty of the MIR spectra between the 
three different epochs. Moreover, because the spectral features are remarkably similar for 
the three spectra in the SL modules, we only focus on the longer wavelengths, from 14.5 to 
32 /im. Finally, due to potential rapid variations with respect to the source's X-ray or radio 
activity, we separately study nod 1 and nod 2, in contrast to using nod- averaged spectra 
for the continuum fitting. All the non-dereddened spectra are displayed in Fig. M and the 
features listed in Table El Note that we only list the features detected at 5a or more above 
the continuum (assessed in 2 fim windows using a second-order polynomial). 

As seen in Fig. [Hland Table |H1 the MIR spectra of Cygnus X-1 are dominated by H I 
and He II emission lines that likely arise from the stellar wind of HD 226868. Overall, more 
features are detected during Obs. 1 and Obs. 2, probably because of a higher soft X-ray 
emission flux level that photoionizes the wind. Moreover, in both cases, the spectroscopic 
content appears to vary from nod to nod, with more lines in nod 1 than in nod 2. This 
is likely due to slight background variations between both positions, although dramatic 
changes of the photoionization state cannot be excluded. 

During Obs. 3, the same features are detected at both positions with the exception of 
two absorption lines at 15.60 ± 0.01 and 25.90 ± 0.02 /im in the nod 1 spectrum, which 
correspond to [Ne III] and He II. A forbidden absorption feature of Ne III is very unlikely, 
and these two lines are probably a consequence of the sky subtraction. Indeed, given 
that slit positions differ between observations, the slit during Obs 3 may be aligned with 
non-homogeneous structures, leading to background variations between the two nods. 

Finally, a total of three blue-shifted H I emission lines are present in the nod 1 and 
nod 2 spectra of Cygnus X-1 during Obs. 2 but not during Obs. 1 nor Obs. 3. Their 
velocities vary between —1600 to —2200 km s~^, which hints at an origin from the stellar 
wind of HD 226868. This is likely due to the orbital configuration of the system during each 



-17- 



observation. Indeed, iHanke et al.l (|2009| ) showed that one should expect the presence of 
blue- or red-shifted X-ray absorption lines from the stellar wind, depending on the orbital 
phase of the BH (Fig. 14 in their paper). Close to phase 0.5, these lines are almost "frozen" 
while they are maximally blue-shifted around phases 0.2—0.3. During Obs. 1, Obs. 2, 
and Obs 3, the BH phases are 0.43, 0. 27, and 0.58, respec tively. Our results are therefore 



consistent with the prediction given in 



Hanke et al. 



(120091). 



6. Discussion 

The main results of our multiwavelength study of Cygnus X-1 are: (1) the compact 
jets contribute significantly to the MIR emission, with a spectral break ^ 2.9 x 10^^ Hz 
(Obs. 1), (2) when the jet carries more power, the cut-off frequency shifts as expected 
to higher frequencies and becomes undetectable in the MIR (Obs. 2), (3) the strength 
of the bremsstrahlung from the stellar winds is likely anti-correlated to the soft X-ray 
emission, which confirms the anti-correlation between the mass-loss and accretion rates (all 
observations) . 



6.1. Detection of the spectral break during Obs. 1 



This, to our knowledge, is the first measurement of the spectral break of a compact 



jet based on spectral fitting instead o 



a break in a BHXB after GX 339-4 flCorbel fc Fender 



photometry, anc 



only the second detection of such 



2OO2I ). We can use the theoretical 



expression of the expected cut-off frequency, to see if our value is consistent with that of 
GX 339—4, scaled with respect to the BH mass, distance, and accretion rate. In the HS, if 
Fi, is the flux density at a given frequency y, M the accretion rate, an d My the BH mass, we 



have for a radiatively inefficient jet, from 



Falcke fc BiermannI (119951 ) and 



Heinz fc Sunyaev 
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fl2003h : 



M2/3 

i/b oc 



X 



(4) 



where i^b is the cut-off frequency. 

For GX 339—4, we consider a distance of 8 kpc, a BH mass of 10 M(^, and a f 



density at 15 GH z of 15 mJy, respectively (jZdziarski et al 



Corbel et al 



2004 



Fender et al 



1997 



ux 



20031 ). For Cygnus X-1 the numbers are 2 kpc, 11 Mq (ICaballero-Nieves et al. 
20091 ) ■ and 15 mJy, respectively. Using Eq HJ we find z/qx ~ 4.06 x i^cyG) where z/qx 
and z/cYG are the GX 339—4 and Cygnus X-1 cut-off frequencies, respectively. For our 
measured vcyg ~ 2.9 x 10^^ Hz, this rela tion gives vcx ~ 1 .15 x 10^^ Hz, which is in 



good agreement with the values derived in 



Corbel fc Fender! (120021 ). Therefore, despite the 



intrinsic uncertainties associated with the MIR spectral fitting, we believe that we see the 
cut-off frequency in the Cygnus X-1 MIR spectrum during Obs. 1. 



6.2. Implications for the jet geometry 



We can us e the presenc e of the synchrotron break to constrain the jet geometry. As 



pointed out in iHeind (|2006| ). the radio photosphere is located at a de-projected distance 
of about 5 X 10^^ cm from the black hole, based on the fact that about 50% of the flux 
are unresolved at an angular resolution of 3 mas in the VLBA observations reported by 



Stirling et al. 



fl200lh . 



For a Blandford-Koenigl jet ( IBlandford &: Koenigl 



19791 ). the location of the 



photosphere along the jet, Zr=i is inversely proportional to the frequency, Zt-=i,v oc v~^. 
Given the marginally higher radio flux in our observation compared to the Stirling et 
al. (2001) observations, the radio photosphere should be located at distance of about 
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6 X 10^^ cm. At a frequency of z/b ~ 2.9 x 10^^ Hz, the IR photosphere would be located at 
a distance of 2;r=i, 2.9x1013 hz ~ 1-7 x 10^° cm. Since the break from optically thick to thin 
occurs when the photosphere reaches the base of the jet, our result would suggest that the 
base of the jet has a scale of about half a light second. For a BH of mass M ^ 10 Mq, this 
corresponds to about 10,000 gravitational radii, which is fairly large compared to size scales 
typically assumed for the base of the jet, usually assumed to occur on scales of a few tens 
to hundreds of gravitational radii. For example, the jet in M87 is resolved down to scales of 
less than 100 gravitational radii. 

A simple Blandford-Koenigl model is a severe oversimplification over the almost 
four orders of magnitude in frequency and thus photospheric scale between IR and radio. 
Deviations are almost certainly to be expected. However, any change would still have to 



obey the observed flat SED. As argued in 



Heinzl (120061 ). simple geometric changes to the 



Blandford-Koenigl model do not alter the behavior of the jet significantly. For example, 
allowing the jet to be continually coUimated (with i?jet oc z^) instead of being conical 
(C = 1) and requiring energy conservation such that p cx 2;"^^ results in a slight altered 
dependence_oj[z3-=Lj^_oc u~^/'^, and a change of the spectral index from completely flat to 



a = 1 — i (jHeinzll2006l ). Bringing the base of the jet inwards (toward the BH) would require 
a C > 1; implying a steeper spectrum, while microquasars tend to have fiat to slightly 
inverted spectra. A simple geometric explanation for the large implied scale for the base of 
the jet is thus not plausible. 

It is clear that the jet must be accelerated between the base of the jet and the radio 
photosphere, changing the Doppler parameter S along the jet. Since the optical depth to 
synchrotron-self absorption is proportional to 5^, a significant increase in 6 from the base of 
the jet, where the IR emission originates, to the radio photosphere, would boost the radio 
optical depth, effectively moving the photosphere outward by a factor of 5^^^. However, 
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given that the jet in Cygnus X-1 is hkely not ultra- relativistic ( IGleissner et aLll2004j ). it is 
unhkely that this could bring the base of the jet inward by more than a factor of a few. 

Finally, it is pos sible that the estimated resolved flux fraction of the radio of ~ 50% 



was overestimated in iHeinzl (j2006l ). For example, if instead of 50%, only 25% of the flux 
were actually resolved in that observation, it would decrease Zr=o by an order of magnitude, 
again moving the size of the IR photosphere and thus the base of the jet towards scales 
more typically associated with jet acceleration. However, very long baseline observations 
tend to underestimate the resolved flux, rather than the unresolved flux. 

If the size scale over which the jet forms is really this large, it would imply that large 
scale outflows (i.e., winds) from the disk are fundamentally involved in jet formation. 
Indeed, a recent study of GRS 1915+105 have s hown that the accre tion disk wind and the 



jets were competing for the same matter supply (INeilsen fc Lee 



20091 ) . Moreover, it has long 



been speculated that winds support jet collimation by providing lateral pressure support. 
This would make the jet in Cygnus X-1 qualitatively different from AGN jets which have 
been imaged down to scales significantly smaller than this in terms of gravitational radii. It 
would also imply that the jet originates over a region much larger than the X-ray emitting 
corona that is typically assumed to be the base of the jet. This underscores the need for 
further confirmation of the nature of the outflow responsible for the observed radio emission 
and, by extension, the large scale diffuse radio and emission line nebula around Cygnus X-1. 



6.3. Origins of the radio and X-ray/7-ray emissions during Obs. 1 

Fig [To] displays the contributions of the stellar photosphere (blue), the bremsstrahlung 
from the expanding stellar wind (magenta), and the compact jet extrapolated to the 
high-energy (green), from our fit of the Cygnus X-l's SED in the case 2, scenario (a). 
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The simultaneous Ryle, Spitzer, and RXTE data are also superimposed, a s well as the 



non-c ontemporaneous INTEGRAL/lBlS data that were recently published (iLaurent et al. 



20111 ). The authors detected for the first time polarized 7-ray emission beyond 400 keV, 
and they argue that it is due to compact jets. Our result is fully consistent with this 
scheme, as the optically thin part of the jet perfectly describes the INTEGRAL/IBIS data, 
providing strong evidence that we actually detect the spectral break around 2.9 x 10^^ Hz. 
Nevertheless, it is clear that the RXTE spectrum of Cygnus X-1 cannot be described by the 
compact jet alone, and one or several extra components are needed. The spectral similarity 
between the spectral indices of the IR and the X-ray spectra might suggest that one of 
these co mponents is synchrotron- s elf-compton (SSC) e mission in the base of the compact 



jet itself (IMarkoff fc Nowak 



2004 



Markoff et al. 



2OO5I ). However, this would require a 



large compactness of that region, in direct contradiction of our discussion regarding the 
implied large size of the IR photosphere and thus the base of the jet above. The more 
natural interpretation would be Comptonization of disk photons by the corona, independent 
of the optically thin synchrotron spectrum from the jet. Moreover, whether or not this 
result implies a disk origin of the X-ray emission in other sources like GX 339—4, whose 
NIR spectra are more consistent with a synchrotron origin of the X-ray, remains an open 
question. 



6.4. A clumpy stellar wind? 



Bremsstrahlung is found to be variable and to contribute significantly to the MIR 
emission of Cygnus X-1, more than the compact jet itself. The fact that its intensity is 
similar in the HS (Obs. 1 and Obs. 3), and lower in the failed-transition state (Obs. 2) 
is consistent with an anti-correlation of the mass-loss and accretion rates, as it was first 
shown by 



Giesetal 



(l2003l ). who derived a mass-loss rate of about 2.57 x 10 ^ Mq yr ^ 
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and 2.00 x lO"*' Mq yi-^ for the HS and SS, respectively. We can also derive from the 
theoretical expression which links it to our fit-derived brernsstrahlung flux dens ity S^, at 



15 GHz, including clumping effects ( iLamers fc Waters 



1984 : 



Scuderi et al. 



19981): 



2.2610"^ X V7^x 



Si/ 



V 



10 GHz 



-0.6 



104 K 



-0.1 



1 kpc 



100 km s"^ 



(5) 



where /oo is the clump filling factor of the wind, To is the electron temperature far from 
the photosphere, which is well approximated as O.ST^,, yU is the mean atomic weight 
of the expanding plasma per electron, usually taken at about 1.3 for O supergiant 



Lamers fc Leit 



GiesetaL 



nerer 



19931 ). and v^o is the terminal velocity of the wind (~ 1600 km s ^ 



20081 ). When clumping is not included (i.e. /oo = 1), we derive mass-loss rates 
of about 8.79 x 10"^ M© yr'^ and 5.95 x 10"^ Mq yr'^ for Obs. 3 (HS) and Obs. 2 (IS), 
respectively. Neglecting clu mping clearly lead s to ov erestimate the mass-loss rate by a 
factor of 3—4, as shown by Lamers fc Waters (1984). This is therefore a hint that the 
Cygnus X-1 w ind is clumped, a nd we need /oo ~ 0.09 — 0.10 in Eq. |5]to match the mass-loss 
rates given in 



Giesetal 



( 120081 ). This value is consistent with that derived from the study 



of the X-ray dips (Hanke et al., in prep). 



7. Conclusion 

We have investigated, through SED fitting, the origin of the MIR emission of 
Cygnus X-1 when the source was in the HS (Obs. 1), the IS (Obs. 2), and in a peculiar 
compact jet-free HS (Obs. 3). We interpret the MIR continuum during Obs. 3 as due to the 
stellar's photosphere and clumpy expanding wind (/oo ~ 0.09 — 0.10). In contrast, when 
jets are present, they clearly contribute to the MIR emission. During Obs. 1, we detect the 
spectral break (z^b ~ 2.9 x 10^'^ Hz), and we argue that this value is consistent with the 
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base of the jet being very large, likely involving accretion disk wind in the formation and 
coUimation of the jet. We show that if the jet alone can account for the 7-ray emission 
of the source, additional processes are needed to explain the 3 — 200 keV spectrum. A 
compact jet is also present during Obs. 2 but we do not detect the cut-off frequency as it 
has likely shifted towards higher frequencies. Nonetheless, the bremsstrahlung emission is 
found lower than that during Obs. 3, which is a confirmation of the anti-correlation between 
the mass-loss and accretion rates. 

Our results highlight the importance of the MIR and, by extension, the far-infrared 
in the study of microquasars, as all their components may have a spectral signature in 
these domains. The existence of powerful infrared space-based observatories, Herschel 
today and the James Webb Space Telescope (JWST) in the near future, combined to 
ground-based telescope such as the Atacama Millimeter/sub- millimeter Array (ALMA), 
are therefore a unique opportunity to make strong progresses in our understanding of their 
behavior and properties. We definitively recommend further studies of Cygnus X-1 and 
other microquasars with these facilities. 
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Table 1. Summary of all the observations of Cygnus X-1 we made use of in this study 

Obs. # IRS RXTE Ryle Ryle fluxes (mjy) $ 

1 53513.0399839-53513.0544390 53513.0433283-53513.2142542 53513.0398560-53513.0546265 5-28 0.426 

2 53528.9592175-53528.9736728 53528.9698098-53529.2433283 53528.9593506-53528.9737549 4-42 0.269 

3 53553.1129871-53553.1274610 53552.9079579-53553.1988838 53553.1131592-53553.1275635 0-48 0.583 



Note. — We give the observation number, the IRS day of observation (in MJD), the day of simultaneous coverage wit h RXTE and the 
Ryle te lescope, the Ryle flux level at 15 GHz (in mJy), as well as the orbital phase 3> calculated from the ephemeris given in lBrocksopp et al. 
(1999a). 
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Fig. 1. — Dereddened spectra of Cygnus X-1 based on Spitzer /IRS spectroscopy taken in 2005 
May 23 (Obs. 1), June 7 (Obs. 2), and July 2 (Obs. 3). 



Table 2. Best parameters from the fits of Obs. 1, Obs. 2, and Obs. 3 i?Xr£;/PCA+HEXTE spectra of Cygnus X-1 



Obs. # ATh Ecut -E-fold Tl Ebreak r2 -Biron firon -^3-20 ^^20-200 X'^/dof 

(10^^ cm-^) (keV) (keV) ••• (keV) ••• (keV) (keV) (10~^erg cm-^ s~^) (10~*erg cm-^ s~^) 



1 O-GSf-iTs 26.15l?°9 175.971^:?^ ^-SStlol W.2it'if6 ^-^Girol 6.381":°^ 0.60l°;«| 1.19 2.41 1.02/303 

2 0.53 (frozen) 26.87t2;Jg 151.67lg;fg 2.08l°:oi lO-SSlo.lg l-70to.oi 6.45l°;°3 0.58t° °5 1.55 2.21 1.08/304 

3 1.04t°;ll 25.28tH° 212.181^:56 1.75l°:°i 10.1ll^|| ^-^St^o.ol ^-^^-oZ 0-39ig;g^ 0.71 1.64 1.09/303 



Note. — The best-fit model is phabs(highecutxbknpower+gaussian), and the errorbars are given at the 90% confidence level. 
"^The fluxed are unabsorbed and exclude that of the iron feature 
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obs. 1 



data and folded model 




Obs. 2 



data and folded model 




Obs. 3 



data and folded model 




Fig. 2.— Obs. 1, Obs. 2, and Obs. 3 RXTE /PCA+REXTE spectra of Cygnus X-1 fitted with the 
model phabs(highecutx bknpower+gaussian) . 
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Fig. 3. — Ryle (15 GHz, bottom) and ASM (1.2—12.0 kev, middle) liglit cm'ves, as well as the 
ASM HR2=C/A hardness ratio (top) of Cygnus X-1 between MJD 53500 and MJD 53580. Our 
three IRS spectroscopic observations are marked in dashed lines. 
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Fig. 4. — Hardness- Intensity Diagram of Cygnus X-1 corresponding to our three observations. The 
X-axis is the 1.2-12.0 kev ASM flux in Crab (1 Crab = 75.5 counts/s), while the Y-axis is the ASM 
HR2=C/A hardness ratio. Clearly, Obs. 3 is the faintest and hardest observation, while Obs. 2 is 
the brightest and softest. 
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Fig. 5. — Obs. 1, Obs. 2, and Obs. 3 Ryle light curves (20 s, blue diamonds) during each IRS 
integration. The black segments represent each module, the height being the median radio flux and 
the length the exposure time. The red line is the median radio flux during the whole observation. 



Table 3. Best parameters from the fits of the Obs. 3 MIR spectrum of Cygnus X-1 



Obs. # T* R^/D^ as Ns xVdof 

(K) (_Rq /kpc) • • • (mjy at 15 GHz) 



3 28000 (frozen) 9.94^^;!^^ O.SSiJj;!]; 0.74t^;?^ 1.00/302 



Note. — The model used is a combination of a black body and a power law. 
The errorbars are given at the 90% confidence level. 



00 



Table 4. Best parameters from the fits of Obs. 1 and Obs. 2 radio/MIR SEDs of Cygnus X-1: case 1 



Obs. # R^/D^ as Ng aj xVdof 

(K) (Hg/kpc) ••• (mjy at 16 GHz) ••• (mjy at 16 GHz) 



1 28000 (frozen) 9.94 (frozen) 0.55 (frozen) QMt^f^ O.Olt|5:^2 15.60t4 34 0-96/301 

2 = = = 0.44lg;^^ 0.06tg;g^ 22.00l|;^^ 0.72/301 



Note. — The model used is a combination of a black body and two power laws. The errorbars are given at 
the 90% confidence level. 



Table 5. Best parameters from the fits of Obs. 1 and Obs. 2 radio/MIR SEDs of Cygnus X-1: case 2, scenario (a) 



Obs. # T* R*/D^ ae Ng aji aja Nj xVdof 

(K) (Hg/kpc) • • • (mjy at 15 GHz) • • • (THz) • • • (mjy at IB GHz) 



1 28000 (frozen) 9.94 (frozen) 0.55 (frozen) 0.75t!^;[J^ "^Trnttrt -0-6 (frozen) 15.40i3;3^ 0.81/300 

2 = = = O-^&to.ol -0-Olto.oi 14.80t};J^ -0.6 (frozen) 21.30t|j^ 0.69/300 



Note. — The model used is a combination of a black body, a power law, and a broken power law. The errorbars are given at the 90% 
confidence level. 



Table 6. Best parameters from the fits of Obs. 1 and Obs. 2 radio/MIR SEDs of Cygnus X-1: case 2, scenario (b) 



Obs. # n R*/D^ aB Nb aji aja Nj xVdof 

(K) (R^/kpc) • ■ • (mjy at 15 GHz) • • • (THz) • • • (mjy at 16 GHz) 



1 28000 (frozen) 9.94 (frozen) 0.55 (frozen) 0.66 (frozen) 0.04t°:oi 29.43^2:29 -0-58toJI I5.20ta:ll 0.85/300 

2 = = = 0.44 (frozen) 0.06tg;g} 42.6ltf j^ l-2lto.36 22.04tgJ^ 0.69/300 



Note. — The model used is a combination of a black body, a power law, and a broken power law. The errorbars are given at the 90% 
confidence level. 



Table 7. Best parameters from the fits of Obs. 1 and Obs. 2 radio/MIR SEDs of Cygnus X-1: case 2, scenario (c) 



Obs. # n R*/D^ aB Nb aji aja Nj xVdof 

(K) (R^/kpc) • ■ • (mjy at 15 GHz) • • • (THz) • • • (mjy at 16 GHz) 



1 28000 (frozen) 9.94 (frozen) 0.55 (frozen) 0.74 (frozen) 0.0lt|5;^i 29.27tJ:^? -l.Uto.tt 15-30l}:4i 0.78/300 

2 = = = = 0.0lt°;°} 15.28t}j^ -0.55t°:°o 21.60t|J° 0.69/300 



Note. — The model used is a combination of a black body, a power law, and a broken power law. The errorbars are given at the 90% 
confidence level. 
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Fig. 6. — Fit of the MIR continuum of Cygnus X-1 during Obs. 3, with a stellar black body and 
an expanding wind bremsstrahlung. 
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Fig. 7.— Obs. 1 and Obs. 2 MIR/radio SEDs of Cygnus X-1: case 1. 
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Fig. 8. — Obs. 1 and Obs. 2 MIR/radio SEDs of Cygnus X-1, case 2: (top) scenarios (a), (middle) 
scenario (b), (bottom) scenario (c). 
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Fig. 9. — LL nod 1 and nod 2 spectra of Cygnus X-1 during Obs. 1 (top), Obs. 2 (middle), and 
Obs. 3 (bottom). 



Table 8. List of all detected features for each LL nodding position 
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Note. — We give the name of the feature, its measured (Af^^) and laboratory (A) wavelengths, its equivalent width W, its full-width at half-length {FW H M) , and its flux. 
^ Blue-shifted lines 
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Fig. 10. — Radio to X-ray SED of Cygnus X-1 during Obs. 1 (case 2a). The compact jet (green), 
bremsstrahlung (magenta), and stehar (blue) contributions are displayed. Simultaneous Ryle, IRS, 
and RXTE data are superimposed. We also add the data from the INTEGRAL/IBIS compton 
mode, obtained by combin ing all the IBIS observations of the source, in the HS, between 2003 and 



2010 ( Laurent et al 



2011 



